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Abstract 

' In our work we extend the ideas of the derivation of the chiral effective theory from 

the lattice QCD [1] to the case of the random lattice regularization of QCD. Such 
procedure allows in principle to find contribution of any order into the chiral effective 
f\| | lagrangian. It is shown that an infinite subseries of the chiral perturbation can be 

summed up into tne Born-Infeld term and the logarithmic correction to them. 
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1 Introduction and motivation 

• r- 1 

At low energies strong interactions can be described by effective chiral lagragians. Such effec- 
tive field approach is very important, both for the physics of the hadron-hadron interaction 
and for the concept of the low-energy baryon state. Therefore, there is a strong motivation 
to derive effective chiral theories from QCD. 

The first essential contribution to this topic was made in [2]. Using the method of large N 
expansion, an effective chiral theory was suggested in the form of a series of chiral invariants. 
The second order theory of such type is the well-known Skyrme model [3] of low-energy baryon 
states, a phenomenological unified theory for mesons and baryons where the baryon is treated 
as a topological soliton of nonlinear chiral fields. Another interesting approach was proposed 
in [4] where the Skyrme model was derived from integration of the chiral anomaly. 

All such methods play a very essential role in particle physics and give appropriate de- 
scription of the behavior of the chiral field at low energy. However, to study the chiral field 
near the confinement boundary one must analyze the whole series of chiral perturbations. 
The Chiral Bag Model [5] provides a very good illustration of this problem. In this model the 
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boundary between chiral fields and color fields is specified by hand. Although this gives quite 
good agreement with experimental data, it is not clear what might be the physical mechanism 
of formation of this chiral bag. 

Fortunately, Lattice QCD numerical experiments give us today a lot of interesting infor- 
mation about the behavior of color fields (quarks and gluons) in the strong coupling regime. 
These data are very essential to understand the physics of the baryon at low energies; in fact 
the lattice approach has formed empirical and theoretical basis of the baryon string model [6]. 

Against the background of these facts, it looks astonishing that the Skyrme model should 
give a good agreement with the experiment data without accounting for the color degrees of 
freedom at all. Therefore the question of unification of these paradigms: the chiral soliton 
model and the model of baryon string, becomes very essential. It is natural to begin the 
search for the way of such unification with analyzing the chiral limit of QCD on the lattice. 

2 Why do we need the Random Lattice QCD 

Derivation of a chiral effective Lagrangian from lattice QCD has been attempted many times 
since long ago. The well-known Brezin & Gross trick [7] makes it possible to perform integra- 
tion of the link matrix in the strong coupling regime and to obtain various first order chiral 
effective theories [8]. 

Although at first such approaches led to great success, they have not been very popular, 
because they do not allow to obtain any corrections to first order results. Lattice regularization 
breaks the rotational symmetry of the initial theory from the continuous rotation group down 
to a discrete group of rotations at fixed angles. Hence, lattice regularization approaches 
give correct results only for those tensors that are invariant with respect to such discrete 
groups. In particular, using the ordinary Hyper-Cubical (HC) lattice, one can obtain only a 
first order effective theory, while for corrections this method generates non-rotation invariant 
(non-Lorentz invariant) terms. Generation of high-order effective field theories requires a 
more symmetrical lattice. 

The problem of breakdown of rotational symmetry on a lattice has been attracting impor- 
tant attention for a long time. It was shown [9] that in 4 dimensions the so-called Body Cen- 
tered Hyper-Cubical (BCHC) or F4 lattice has the largest discrete symmetry group. (BCHC 
consists from the all sites of the HC lattice together with centers of its elementary cells.) This 
property of the BCHC lattice gives a possibility to obtain the next-to-leading (NL) correction 
to the first order of the chiral perturbation theory [1]. 

The results of the papers [1] are essential for our analysis, as they confirm the effectiveness 
of the idea of chiral effective lagrangian derivation from the lattice QCD. Moreover, these 
results are interesting from phenomenological point of view because, as is well known [3], 
the NL corrections violate the scale invariance of the prototype (first order) chiral theory 
that leads to generation of chiral topological solitons (Skyrmions). The Next-Leading order 
chiral effective theory that was implemented in [1] is in agreement with our phenomenological 
propositions [10], and in our work we will use methodological ideas from [1] in order to define 
the behaviour of chiral field near the confinement surface. 

As one could see, in order to solve our problem, the Next-Leading order corrections are 
not enough. This theory has no solutions that look like chiral "bag". Moreover, as will be 
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shown later, near the confinement surface (near the source of the chiral field), the influence of 
high order corrections became larger and larger. But the BCHC lattice method gives the NL 
corrections only and the further use of this method for the defining of the high order terms 
leads to generation of non-relativistic (non-rotational) invariant terms. It means that we need 
a more symmetric lattice than the BCHC lattice. 

Unfortunately, a lattice which would be more symmetric than BCHC lattice cannot be 
constructed in 4 dimension. Moreover any method based on a lattice of a fixed geometry 
has artifacts coming from priority directions that correspond to basis vectors of the lattice. 
It is these artifacts that eventually lead to the problems with the rotational (relativistic) 
invariance rendering the use of the BCHC lattice to be only half measure. For solving our 
problem a modification of the initial concept of lattice regularization must be performed. We 
need to find a concept of lattice regularization that has no priority directions. Fortunately 
this concept is known for a long time and is called the Random Lattice approach [12]. 

The idea of Random Lattice was proposed originally by Voronoi and Delaunay: today this 
method is widely used in the modern science. For the quantum field theory the method was 
modified by Christ, Friedberg and Lee [12]. In these articles it has been shown that in order 
to restor the Lorentz (rotational) invariance, it is necessary to perform an average over an 
ensemble of random lattices. As a result one gets the averaging over all possible directions 
and it is intuitively clear that this procedure leads to the disappearance of the artifacts that 
cause the violation of the group of the space rotations. 

But how to perform such random discretization? This procedure has the tree steps: 

1) Pick iV sites Xi at random in the volume V. 

2) Associate with each Xj a so-called Voronoi cell q 

Ci = {x | d(x,Xi) < d(x,Xj),Vj 7^ i} 

where d(x, y) is a distance between points x and y. It means that the Voronoi cell Q consists 
of all points x that are closer to the center site x^ than to any other site. 

3) Constrict the dual Delaunay lattice by linking the center sites of all Voronoi cells which 
share a common face. 

Now if one considers the the big ensemble of such Voronoi-Delaunay random lattices based 
on various distributions of sites Xi, it possible to prove that the original rotational symmetry 
is restored [12]. In our work we use this procedure to obtain an effective chiral lagrangian from 
lattice QCD. This methodological point of view it is a modification of the method proposed 
in [1] in the case of the Random Lattice approach. 

3 From Lattice QCD to chiral lagrangians: step by step 

Now let me briefly recall a general steps of the algorithm of derivation of the chiral lagrangian 
from the lattice QCD that was proposed in [1]. 
Step 1: Definitions 

The starting point of our analysis is a standard lattice action with Willson fermions 
Z = J [DG] [Dip] [Dip] exp{-S pl (G) - S q (G,$,i/,) - Sj} 
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where: 

1) the plaquette gauge field term is defined by 



2N C r 1 i r 

pi 



2)the link fermions term is defined by 

S q = Y,to(Mx)Gv{x) + G+(x)Mx)) 



A^{x) a b = Mx + /i)P+^ a (x), A,{x) a b = i, b {x)P~r{x + fj) 

and P± = |(r ±7 M ); 

3) the source term is defined by 



x JV c 



Step 2: Strong-coupling regime on the lattice and integration over the gauge 
field 

In order to realize the strong-coupling regime on the lattice let us consider the limit of the 
large coupling constant g (g — > oo). This limit was widely studied [8] and the main result is 
that in such limit integration over the gauge field can be performed. (Of course, the direct 
integration is difficult since there exists the plaquette term S^, but due to the strong-coupling 
limit on the first step plaquette contributions are negligible with respect to the contribution 
from the link integral S q . The plaquette contributions could be considered in the systematic 
manner as perturbations in 1/g [8].) 

Let us consider the leading order contribution in this strong-coupling expansion. The 
integrals over the gauge degrees of freedom can be calculated into the large N limit by using 
the standard procedure [8] and the result of these calculations is the following 

Z= [[D^[D^}exp{-Nj2t4F(Hx,v))}-Sj}, (1) 

J x,v 

where A„ = -M(x)P~M(x + u)P+ and 



F(X) = tr[(l - v / T^A)] - tr[log(l - -VT^X)}. 

Interestingly, the function F(X) has the typical form of the Born-Infeld action with a first 
logarithmic correction. This is no coincidence. In [11], it was shown by means of very similar 
technique that the low-energy theory of the IIB superstring has a Born-Infeld form. From the 
methodological point of view we perform a similar analysis for QCD on the lattice and it is 
important to note before starting our proof that our result will have a Born-Infeld form too. 

Step 3: Integration over the fermion field and chiral limit 

Our next step is the integration over the fermion degrees of freedom in (1). Using the 
source technique it was shown [1] that integral (1) can be re- written into the form of an 
integral over the unitary boson matrix M x 

Z = J DMexpS eS (M). (2) 
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As a matter of principle, we already performed the transformation from the color lattice 
degrees of freedom (G and ip) to the boson lattice degrees of freedom (M). Now our task is 
to realize the continuum limit of expression (2). 

This step of our analysis amounts to studying of the stationary points of the lattice action 
^eff" Fortunately this is a very well studied task [13]. This problem is connected with well- 
known investigations of the critical behavior of the chiral field on the lattice and with the 
problem of the phase transformation on the lattice (for references see issue [14]). In [1], it 
was shown that for our task the stationary point is 

M = u I, u (m q = 0,r = 1) = 1/4. 

Now one can expressed M(x) in terms of the pseudoscalar Goldstone bosons 

1 ~t~ T5 1 — T5 

M = w exp(i7r i r i 7 5 // 7r ) = u [U(x) — h U + (x) — - — ] 

and the effective action is given in the form of the Taylor expansion around this stationary 
point 

S eS (U) = -N± ^^Etr^Or) - A ) fc ]. (3) 

fc=l x,u 

Let us consider the expansion of the chiral field U = exp(i7rjr i // 7r ) on the lattice around 
the point x in power of the small step of the lattice a 

a 2 

U(x + v) = U(x) + a{d u U{x)) + —(dlU{x)) + ■■■. 

And for components of the Taylor expansion (3) one obtain 

tr[(A„(a;) - A )] = -2A tr(a) 

ti[(\„(x) - A ) 2 ] = 2A 2 tr(a 2 ) -4A 2 tr(a) 

tr[(A,(x) - A ) 3 ] = -2A 3 tr(a 3 ) +6Aj]tr(a 2 ) , (4) 

tr[(A„(a;) - A ) 4 ] = 2A^tr(a 4 ) -8A 4 tr(a 3 ) +4A 4 tr(a 2 ) 



tr[(X u (x) - A ) 5 ] = -2A^tr(a 5 ) 

where a = a?d v Ud v U + + 0(a 4 ). 

Step 4: Problem of rotational symmetry violation: examples of the Hyper- 
Cubical and Body Centered Hyper-Cubical lattices 

Expressions (4) are very essential because these are a simplest illustration of all aspects of 
the violation of the rotational symmetry on the lattice. For this moment we assume nothing 
special about the structure of our lattice. We try to formulate our result as generally as 
possible and all information about the lattice contained in the vectors v that correspond to 
the basic vectors of the lattice (for example, the vectors v for the Hyper-Cubical lattice are 
the Cartesian basic vectors i, j, k and t). The leading order part can be calculated trivially. 
Indeed, using the simple Hyper-Cubical lattice where v — i,j : % — (1 . . . 4) it is easy to show 
that the leading order contribution is the prototype chiral lagrangian 

Po(p 2 )~tr[^C/^C/ + ]. (5) 
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As I said before the rotational symmetry violation argument does not allow to use the HC 
lattice calculation for the next-leading order contributions. For obtaining of these contribu- 
tions a more symmetrical lattice must be used. In [9] it was shown that this lattice is a Body 
Centered Hyper-Cubical (BCHC) or F4 lattice. 

Unfortunately, this method can not be directly used for finding next contributions and 
the origin of this fact is again the violation of the rotation symmetry but now on the F4 
lattice. Moreover, there are no any more symmetrical lattice with fixed positions of sites in 4- 
dimensional [9]. It means that we need an absolutely different lattice concept that guarantees 
the restoration of the initial symmetries. Fortunately this concept is known now. This is the 
Random Lattice concept (RL) [12]. 

4 Random Lattice in action 

The basic idea of the RL is the averaging over the big ensemble of various lattices with random 
distributions of sites and it is possible to show that such averaging leads to the restoration of 
the rotational invariance. There are two methods of the realization of such scheme. A first 
one based on the Christ, Friedberg and Lee (CFL) technique [12]. 

Commonly CFL technique leads into complicated geometrical analysis. For our task it 
would be very useful to use the analogy between Random Lattice and Random Surface tech- 
nique that was revealed recently [15, 17]. The idea is quit simple: for beginning let us consider 
a lattice with fixed positions (for simplicity it is possible to use the trivial HC lattice, where 

— * — * 

basis vectors are just V — i,j . . .) in a flat space. For a simulation of the Random Lattice 
let us consider small deformations of the geometry of this space (7^ — > gij) so that one can 
rewrite the problem of the random lattice averaging in the terms of the random deformations 
of the geometry of this space [15]. This is a standard quantum gravity problem for which 
powerful methods of the Matrix Theory could be used. 

In our problem we discuss the link integrals that depend on the basis vectors v. All such 
integrals are considered separately for any lattice site X{. It means that rotation invariance 
violation artifacts could be avoided by considering only rotation deformations of these basis 
vectors (translation and re-scaling deformation are left aside in our case). 

Let R G SO (4) be a rotation operator of 4 dimensional vectors v 



It is essential to note that our task can be reformulated in the language of the standard 
Hermitian averaging because SO(4) = SU(2) x SU(2). For infinitesimal rotations we obviously 
have 



where is a traceless antisymmetric matrix. 

Consider the Gaussian Ensemble of such arbitrary rotation. The matrix average of an 
arbitrary function / with respect to Gaussian measure is 



v\ - RijVj . 




(6) 
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where dH is the standard Haar measure. The normalization factor Nq is fixed by requiring 
that (/ = 1)3 = 1. 

Using of the Matrix integration technique one can prove the so-called Matrix Wick's 
theorem for traceless antisymmetric matrices [18, 15, 16, 17] 

{HijH k i)^ = SijSji — SuSjk, 

(U h v)s!>= e (-irnw^. ( 7 ) 

(«>j) pairings pairs 

where the sum extends over all possible pairings and k is the number of crossings in the 
pairing. The matrix average of any odd combination of ify equals zero due to the parity 
argument. 

It is not hard to prove that the main contribution into the averaging sum over such 
infinitesimal rotations comes from the original non-deformable lattice (connected with the 5ij 
part in (6)). In order to cancel non-deformable lattice artifacts let us consider the ensemble 
fj without this non-deformable contribution: Sj' — ft — 0. Using such averaging, for leading 
order of chiral effective lagrangian one gets 

(tr(a))^* = (tY{a 2 d v Ud v U + ))fi = tr(3a 2 d i Ud i U + ) = tr^a 2 ^^)), 

where Lj = U + diU, and for the NL correction one obtains 

(tr(« 2 % = (ti(a A d v Ud v U + d v Ud v U + ))tf = tr(3V((L,L,) 2 - ^[L t , L 3 ] 2 )) = . 

= -3V Qtr 2 (LjLj) + tr 2 (LiLj) - 4tr(L i L i L J L j )) = (8) 

These results reproduce the HL and BCHL results and it is easy to see that this is just what 
we expected to receive because this contribution was obtained from many other approaches 
[10]. In such a way, we can apply this procedure to corrections of any orders from (4) and 
obtain the rotation invariant result due to the pairing. It means that the question about the 
derivation of the chiral effective lagrangian from Lattice QCD become just a combinatorial 
task. 

It is interesting to point out that our SU (2)-flavor result for coefficients in the NL order 
contribution of the Chiral Perturbation Theory (8) 

Lr t r 

lj 2 _ -^3 

1 2 4 

is in agreement with experiment data from ttti — > titt scattering [10] and these coefficients 
torn out to be closed to the model prediction of the V exchange [19]. 

In the last part of this section we show an application of our procedure. We will find that 
an infinite subseries of the chiral perturbation can be summed up into the Born-Infeld form. 

If the expression (7) allows us to calculate all terms in expansion (4), let us consider just 
the first column there. It is easy to show that either of these is proportional to some power 
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of the leading order contribution (5) 

tr[a] 

tr[a 2 ] r 
tr[a 3 ] r 



tr[(L^) 3 ] 



+ 
+ 
+ 



(9) 



tr[a n ] ~ ti[(L^) n ] + 



Substituting (9) into (3) and collecting all terms which depend on the power of the prototype 
lagrangian one obtains the following expression for the effective chiral lagrangian 



-eff 



-tr 



1 - yfl - - tr [log(l - hi - y/l - 1/(3 2 L^)) 



+ 



(10) 



where • • • are all other terms (in particular the Skyrme term) and j3 is an effective coupling 
constant that depends on the value of our stationary point uq. 

Now let us discuss the result (10). It was obtained that some part of chiral effective action 
has a Born-Infeld form plus a first logarithmic correction to it. In [20], it was shown that 
such form of the effective action plays a very essential role in the problem of the chiral bag 
formation because just these square-root terms generate the step-like distribution solutions 
that can be interpreted as internal phases in the two-phase model of the low-energy baryon 
states. Another terms play essential role only on large distances from the confinement surface 
and can be considered as corrections. 



5 Conclusions 

The aim of this paper is to derive the chiral effective lagrangian from QCD on the lattice 
at the strong coupling limit. We find that this theory looks like a Born-Infeld theory for 
the prototype chiral lagrangian. Such form of the effective lagrangian is expected. From the 
methodological point of view our consideration is very similar to the low-energy theorem in 
string theory that leads to the Born-Infeld action [11]. Moreover, in [20], it was shown that 
Chiral Born-Infeld Theory (without logarithmic corrections) has very interesting "bag" -like 
solutions for chiral fields. It was an additional motivation of our work. 

The Chiral Born-Infeld theory is a good candidate for the role of the effective chiral theory 
and a model for a chiral cloud of baryons. In this model one can find not only spherical "bags" , 
it is possible also to study the "string" -like, toroidal or "Y-Sign"-like solutions, or some other 
geometry. The geometry of the confinement surface depends directly on the model of color 
confinement and it would be very interesting to use, for example, the Lattice QCD simulations 
for the color degrees of freedom in combination with our model for the external chiral field. 
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